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Abstract 

Whitham theory of modulations is developed for periodic waves described by non- 
linear wave equations integrable by the inverse scattering transform method associated 
with 2x2 matrix or second order scalar spectral problems. The theory is illustrated 
by derivation of the Whitham equations for perturbed Korteweg-de Vries equation and 
nonlinear Schrodinger equation with linear damping. 

1 Introduction 

The theory of perturbations of soliton solutions of integrable equations is well developed 
and has found many applications to various physical problems (see, e.g. ^). On the con- 
trary, perturbations of wave trains did not attract much attention, although there are many 
problems which should be treated in framework of such perturbational approach. As the 
most evident applications we can mention the theory of dissipationless shock waves in non- 
linear media described by the Korteweg-de Vries (KdV) [2], nonlinear Schrodinger (NLS) 
jH], derivative NLS ji], Kaup-Boussinesq (KB) jH], and many other equations. As other 
examples one can also mention fluxon trains in long Josephson junctions 6J described by 
the sine-Gordon equation and soliton trains in Bose-Einstein condensate confined by mag- 
netic traps [7] and described by the NLS equation with external potential. In all these 
cases unperturbed nonlinear equations have periodic solutions and typical problems can be 
often reduced to finding slow evolution of parameters of the nonlinear wave due to its initial 
modulation and small perturbation terms in the governing wave equations. 

The general approach to treatment of such problems was suggested long ago by Whitham 
jHJ in] and it is based on two scales of time — fast oscillations of the wave and slow evolution 
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of its parameters. Then equations for slow evolution of the parameters can be averaged 
over fast oscillations what leads to the system of the first order partial differential equations 
called Whitham equations. Whitham discovered that in the case of (unperturbed) KdV 
equation the arising in this way system can be drastically simplified by introduction of 
proper parameterization of the periodic solution when Whitham equations take the so-called 
diagonal form for slowly varying parameters called Riemann invariants. Later it was realized 
pn] that this property of the KdV equation is connected with its complete integrability by 
the inverse scattering transform (1ST) method and just this method yields in a natural 
way the periodic solutions of the KdV equation parameterized by the Riemann invariants. 
In the fundamental paper ^Uj Whitham equations were derived for the general case of 
quasi-periodic solutions of the KdV equation. After that this method was applied to other 
integrable equations El E] and used for description of dissipationless shock waves in 
modulationally stable systems jEl 1^ IHj and formation of soliton trains at sharp fronts of 
pulses in modulationally unstable systems (see ^H] and references therein). 

The Whitham theory approach implies that the process under consideration is considered 
at asymptotically large periods of time much greater than the period of wave oscillations. 
However, it is clear that small perturbations neglected in completely integrable approxima- 
tions can make considerable contribution into long-time evolution. Therefore it is necessary 
to develop perturbation theory in the Whitham approach and the first step should be devel- 
opment of the method of derivation of perturbed Whitham equations. The method suggested 
in [T7j for perturbed KdV equation was quite complicated and has not found many appli- 
cations. Therefore the authors of refs. [T8| IT9| 120] used the classical Whitham approach 
for consideration of influence of weak dissipation on long-time evolution of the KdV wave 
trains. In a similar way, the Whitham equations for two-phase quasi-periodic solutions of 
perturbed sine-Gordon equation obtained in [21] did not have effective enough form and 
were considered in a very special case only. 

The aim of this paper is to develop simple and general enough method of obtaining 
Whitham equations for modulated wave train solutions of perturbed nonlinear equations. 
The method can be applied to integrable nonlinear wave equations associated, in framework 
of the 1ST method, with second order scalar spectral problem or equivalent to it 2 x 2 matrix 
problem (Ablowitz-Kaup-Newell-Segur scheme [221 and its generalizations). The general 
formulae obtained in the next section are then illustrated by applications to the perturbed 
KdV equation (where known results are reproduced) and to the NLS equation with weak 
linear damping. These examples show that our approach is quite effective and applicable to 
many nonlinear wave equations of physical importance with various types of perturbations. 

2 Perturbed Whitham equations 

Many physically important equations can be solved, in principle, analytically by the 1ST 
method in which they are represented as compatibility conditions of two linear systems with 
spectral parameter (see, e.g. [22] )■ We shall confine ourselves to the cases when these two 
systems can be represented in one of two equivalent forms. The 2x2 matrix form is given 
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by 

where the matrix elements depend on as field variables Um{x,t) of the equations under 
consideration, so the spectral parameter A. For example, in the case of the NLS equation, 

iwt + Mxu' ± 2|m|\ = 0, (2) 



these matrix elements correspond to the Zakharov-Shabat spectral problem |23] and are 
equal to 

F = — iA, G = iu, H = ±i|-up; 

A = -2iA2 ±i|M|2, B = 2iu\-u.^, C = ±2iu*\±ul, 

and there are two field variables, u{x, t) and its complex conjugate t). The compatibility 
condition "^^t = "^tx of the systems ((T)) yields at once the evolution equations in general form 

Ft- A^ + CG - BH = 0, 

Gt-B^ + 2{BF - AG) = 0, (4) 
Ht-C^ + 2{AH - CF) = 0. 

Taking as an example again the NLS equation, it is easy to see that if we substitute © into 
(01) and demand that these equations must fulfil for any value of the spectral parameter A, 
then we find the NLS equation and its complex conjugate, 

Ut = iuxx =F 2i|'up-u, ul = —iul^ ± 2i|Mp'u*. (5) 

In general case we shall write symbolically the equations in terms of the field variables 

t) as 

Um,t = Km{Un,Un,x, ■ ■ ■), TU, H = 1 , . . . , N , (6) 

Many physically important nonlinear equations are associated with linear systems IQ. 

Another common form of linear equations associated with nonlinear integrable equations 
can be written as 

ipxx = Alp, (7) 

iPt = ~B,ij + Bij„ (8) 

where again the coefficients A and B depend on the field variables and the spectral parameter 
A. For well-known example of the KdV equation 

Ut + 6uUx + Uxxx = (9) 

there is only one field variable u{x,t) and 

A=-{u + \), B = A\-2u. (10) 
The compatibility condition {jpxx)t = {4't)xx of the equations 0,® yields 

At - 2B,A - BAx + Ib^xx = 0. (11) 



After substitution of (fTUj) this equation yields the KdV equation and in the general case 
O,® the resulting evolution equations can be represented again in the form (jH)). Besides the 
KdV equation, the representation ((Tj),® is convenient for treatment of the KB system ^25| 
and some other equations. As was shown in [2^] , the matrix system (Q) can be transformed 
to scalar form (j?!),® by the formulae 

A={F- G^llGf + GH + {F- (12) 

B = B/G, (13) 

so that equations (jH) reduce to (fTT|) . Thus, these two forms of linear spectral problems are 
equivalent to each other and both can be used for treatment of the respective nonlinear 
equations. The scalar form (|7j),(jHl) is more convenient for derivation of perturbed Whitham 
equations and therefore it will be used in what follows. 

We assume that unperturbed wave equations (0) have periodic solutions which can be 
obtained in framework of the 1ST method in the following way (see, e.g. [23j). We take two 
basis solutions and ip' of the second order differential equation ((Zj) and construct from 
them the so-called 'squared basis function' 

g = ij+^-. (14) 

It is easy to show that it satisfies the equation 

g^xx - 2Axg - 4:Agx = 0, (15) 

which upon multiplication hj g/2 can be integrated once to give 

lgg..-\gl-Ag' = PiX), (16) 

where the integration constant denoted by -P(A) can depend on the spectral parameter A. 
In the finite-gap integration method (see, e.g. |E2j) it is shown that periodic solutions are 
distinguished by the condition that -P(A) be a polynomial in A . Then as a function of A 
should also be a polynomial in A with coefficients expressed in terms of Um{x, t), and, hence, 
Um{x, t) can be found from equations obtained by equating coefficients of A'^ in Eq. (fT^ (we 
suppose here that ^ is a rational function of A in accordance with majority of applications). 
The dependence of g, and hence of its coefficients, on time t is obtained from the equation 



gt = Bg^ - B^g (17) 

which follows from Eqs. (jHI) and (fT^. As a result, Um{x,t) are expressed as (quasi)periodic 
functions of x and t parameterized by the zeroes Aj of the polynomial -P(A) = nfii(A — Aj), 

Um{x,t) = Um{x,t;Xi,. . . ,Xm), m = l,...,N. (18) 



In strictly periodic solutions the parameters Aj are constant by definition. But we can 
slightly modulate this periodic wave by making Aj slow functions of x and t which change 
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little in one wavelength and one period. Whitham showed jHl El (for particular case of one- 
phase periodic solution of the KdV equation when this solution is parameterized by three 
zeroes Ai, A2, A3) that substitution of 

Um{x,t) = Um{x,t]Xi{x,t),. . . ,\M{x,t)), m = l,...,N, (19) 

into conservation laws of Eqs. (jH)) followed by averaging of these conservation laws over fast 
oscillations of the wave yields for slow variables Xi{x,t) the so-called Whitham modulation 
equations 

^ + t'.(A,)|^=0, t,j = l,...,M, (20) 

which govern slow evolution of modulated periodic solutions of unperturbed equations (0). 
This approach is well developed and finds a number of applications (see, e.g. |2S])- 
Now, we want to take into account perturbations in the wave equations (jUj). 

Um,t = Kjn{Un,Un,x,- ■ ■) + Rm{x,t,Un,Un,x,- ■ m, n = 1, . . . , iV, (21) 

where the perturbation terms Rm can be slow functions of x and t and can also depend 
on the field variables Un and their space derivatives. The perturbation terms can have any 
nature, e.g., they can arise due to external potentials or small physical effects neglected in 
derivation of unperturbed equations (0). Obviously, the functions (fTHj) are not the solutions 
of the perturbed equations (PT|) anymore, but again, as in the case of slow modulation, we 
can describe effects of perturbations as slow change of the parameters Aj, i = 1, . . . , M. Then 
after averaging over fast oscillations of the wave we have to obtain additional terms in the 
Whitham equations (j2(jp which represent effects of perturbations. Our aim is to derive such 
generalization of the Whitham equations. 

As was indicated by Whitham 11113; most convenient way of obtaining the modu- 
lation equations is to average the conservation laws. Completely integrable equations have 
an infinite sequence of conservation laws and averaging of their generating function yields 
for periodic solutions of equations (0 associated with the spectral problems ((H) or (jZj),® 
the modulation Whitham equations directly in diagonal Riemann form ^U]. Therefore the 
starting point of our analysis should be the formula for generating function of conservation 
laws for equations (fTTj) described by the 1ST scheme ((Tj),®. As was mentioned above, we 
suppose that ^ is a rational function of A. Then it can be chosen so that 

^Ia^oo (22) 

where a = ±1 and r is an integer number. With such a normalization of A, we can rewrite 
Eq. (jl6|) in the form 

lggxx-\9l-Ag'' = -aPiX), (23) 
where P and g are polynomials in A, 

M degg 

P = J](A-A,), g=l[{X-fi,{x,t)), (24) 

i=l 1=1 
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such that their degrees are connected by the relation M = r + 2 degg. If we introduce new 
function 

Xr/2 

then the identity ()23p takes the form 

^99xx - -^gl - Af = -a, (26) 

where 

A = ^/A^ (27) 

and 

^Ia^oo = 1- (28) 
Thus, we can expand g in inverse powers of A, 

oo 

^ = Ef' ^0 = 1, (29) 

and the coefficients gn can be calculated in recurrent way by substitution of into (j2EI)- 
The examples of this procedure will be given below. 

Now, from (fTTj) and we find that g satisfies the relation 



(30) 



which can be considered as a generating function of conservation laws. Indeed, the coefficients 
of the expansion 

-. oo 

^ n=0 

can be expressed in terms of the already known coefficients with the use of trivial identity 
g-{l/g) = 1, and then these coefficients g-n serve as the densities of the conservation laws and 
the coefficients of similar expansion of B/g in inverse powers of A serve as the corresponding 
fluxes, 

{g.r.\ = [{B/~g).nl, (32) 
where {B/g)_n denotes the coefficients in the expansion 

B_^{B/gU 

^ n=0 

The functions g and 1/g satisfy important identities 
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derived in the Appendix. Here S/Sum denotes the "Euler derivative" 

6 d d d d'^ d -,\Nmn ^^""^ 



6Um dUm dx dUm,x dx^ dUm,xx ' ' ' dUm"""'' 

when it is apphed to g-n-, Nmn is the order of the highest derivative Um"""^ in g^n, and Am 
is the order of the highest derivative ul^"-^ in A. 

Now let us turn to perturbed equations and calculate the correction to the conser- 
vation law ()32|) with the density 



g—n — g—n{Urm '^Tn,xi ■ ■ ■ i '""''). (36) 

Evidently, we have 

{9-n)t ~ '^^m,k 9-n,mkdtUm^ = "^^jk 9—n,mk{dx)^ dfUm 

^^m,k 9—n,mk{dx) ~\~ ^^^rn,k 9~n,mkipx^ Rmi 

(k) 

where g-n,mk denotes the derivative of g-n with respect to Um . The terms with must 
collect into divergence of the flux, dx[{B/g)-n\, so that we obtain 

{g-n)t = [il3/~g)-nl + Yl 9~nMdx)'Rm. (37) 

m,k 

In what follows we shall average these conservation laws over some interval A of x which is 
supposed to be much greater than the wavelength L and much less than the characteristic 
interval along which the parameters of the wave change considerably. Therefore we can 
transform the terms with Rm by "integration by parts" and neglect the contribution of 
boundaries because in the averaged form it is divided by A and goes to zero in asymptotic 
limit A ^ L of very slow modulation. As a result of this transformation we obtain the 
conservation law in the following equivalent form: 



{9-n) 



After dividing this formula by A" and summation over n we arrive at the generating function 
of perturbed conservation laws 



For further transformation we apply the identity ()34j) and get 




(38) 



Here we again can "integrate by parts" and neglect the boundary terms to obtain 
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To average this expression, it is convenient to return to the squared basis function g with 
the help of Eq. (fTTfl . 



OA d'-Rrr, q ^^Q^ 



Here the zeroes Aj of the polynomial -P(A) are slow functions of x and t which have to be 
differentiated with respect to x and t in the left hand side of Eq. (|4m) . Then we obtain the 
terms with 

1 ^ 1 dX, 

^r— and — — , 
V A - Xi ot VA - Xi ox 

which are singular at A Ai, and in the right hand side of Eq. ()4()|1 we already have terms 
with the same singular behaviour. Since the perturbed conservation laws ()40|) after averaging 
should be true at any A including A,, we obtain 

dXi dXi a / dA \ ■ 1 ,1^ ^^^^ 



where 



^^/^^ . = 1,...,M, (42) 



angle brackets denote averaging over proper interval of x, and everywhere A is put equal to 
Aj. Formula ()42j) gives well-known expressions for characteristic velocities in unperturbed 
Whitham equations. The right hand side of (PT|) describes the effects of perturbations. Equa- 
tions 1)411) represent the result we wanted to obtain. The averaged values of the expressions in 
angle brackets can be evaluated with the use of formulae for periodic solutions of particular 
integrable equations under consideration (see, e.g., [TUl 

Let us illustrate the obtained expressions for perturbed Whitham equations with two 
examples. 



3 Examples 

We shall consider here simple examples of application of general formulae (jl^,(jl21) to one- 
phase periodic solutions of the KdV and NLS equations. 

3.1 Korteweg-de Vries equation 

In the case of the KdV equation 

ut + 6uu^ + u^^^ = R (43) 
we have expressions (fTIH) for A and B. From we find the expansion 

. ^ u/2 1 /3 2 1 \ 1 / 5 3 5 2 5 1 ,jys\ 

. = 1 - ^ + ^ (s" + 8^- ) - A^ ( 16^ + 32^^' + 16^^- + 32^' ^ J + " " " ' ^''^ 
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hence 



- = 1 + TT7 -M + -Uxx + —u H K H uu^^ H u^^^> + . . . , (45) 

g A A2 V8 8 J A3 Vl6 32 ^ 16 32 / ' v ; 

and it is easy to check that these first terms of the expansions satisfy the identity i.e. 

? 1 g 



The conservation laws following from Eq. (jH^ have well-known form (see, e.g. 

The periodic solution ("cnoidal wave") of unperturbed KdV equation has the form (see, 

e.g. m) 

u{x,t) = As - Ai - A2 - 2(A3 - A2)sn2(v/A3-Ai(x - Vt),m), (47) 

where 

V = -2(Ai + A2 + A3), m = ^^^, (48) 

A3 — Ai 

and the parameters Ai < A2 < A3 are the zeroes of the third degree polynomial 

3 

P(A) = J](A - A,) = A^ - s,X' + S2X - ss (49) 

i=l 

defining the solution according to Eq. (fT^ . Then 

g{x,t) = X- ij,{x,t), (50) 
where fJ,{x,t) is connected with u{x,t) by the relation 

M = 2/i-(Ai + A2 + A3) =2/i-si (51) 
and depends only on the phase 



e = x-Vt, /ie = 2v/-P(/i), A2</i<A3. (52) 
The wavelength of the solution (jTTj) is given by the expression 

' U r ^■^('"> (53) 

K(m) being the complete elliptic integral of the first kind. In this one-phase case we can 
average over wavelength and the Whitham velocities are easily expressed in terms of L by 
means of the formula 



1 \ / 1 \ 1 dx If d/i 2 dL 



- ) 



.g/\=\ \Xi-fi/ L Jo Xi-fi 2L J (Aj - /i)A/-P(/i) LdXi' 
so that we obtain 



(54) 
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Substitution of ()53p into this expression yields the well-known formulae obtained by Whitham 
jHl El • In the KdV equation case the sum in the right hand side of Eq. pij) reduces to a single 
term with dA/du = — 1. Taking into account that a = —1, g = \i — n = (2Aj — Si — m)/2, 
we arrive at 

dt ^""'dx- dL/dK 4 n,v.(A. - A,) ' ^ - 1' ^' ^^^^ 

This is a compact expression for perturbed Whitham equations obtained in pUj . 

Further calculation can be done after specification of the perturbation R. For example, 
if we take 

R = eu^^, (57) 

[e is a small numeric parameter) what corresponds to taking into account small viscosity in 
hydrodynamic applications of the KdV equation, then we have 

dP 

u^^ = 2fi^^ = -2—, 2\i- si-u = 2{Xi- fi), 
d/i 



so that 



and standard calculation of this elliptic integral yields 

((2A. - s, - u)eu^^) = [(1 - ^ + m')E{m) ^^^^ 

-(1 -m)(l -m/2)K(m)], 

E(m) being the complete elliptic integral of the second kind. Substitution of into 
()56|) gives the perturbed Whitham equations for the KdV-Burgers equation found first in 
[IHlIini- Thus, the general formula pT|l reproduces quite easily the already known results 
for the perturbed KdV equation. 



3.2 Nonlinear Schrodinger equation 

As the second example, let us consider the focusing NLS equation with linear damping 

iut + Uxx + '^\u\'^u = —^£U- (59) 

(Similar results for the defocusing NLS equation will be presented in the end of this Section.) 
The unperturbed equation is usually associated in framework of the 1ST method with the 
2x2 Zakharov-Shabat spectral problem [21]. However, with the use of Eqs. f|T^ .(fT ^ it can 
be transformed into scalar form (j7l8|) with (see [23120]) 

Expansions of g and 1/g are given by 

, 1 iwx 1 /I^P Uxx\ 1 fiuxxx 3mxU*\ 

' = ' + -X2^-Y^[— + J^)-Y^[-^ 
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1 

9 



1 iux 1 



u 



2~ 4^ 



1 / iWa-xx iMa;W* 

l~8^^ 4 



iu 4 

and it is easy to check that they satisfy the identities 

5 1 _ _ 1 
~5u^ ~ ~2A2 



+ 



ml 



Ux 

4u2 



+ ... 



(62) 



'dA^ d / 


dA\ 






du^ dx \ 


dujj 


dx"^ 


\du^x ) _ 


S 1 


u 






Su* g 









The generating function ()30|) yields the conservation laws in unusual form 



/ i^x 



2\uV + 



u 



2 , ^^xx 



m + 



2m 2m2 



2m2 



2m 



2mt.m* 



(63) 
(64) 

(65) 
(66) 



The one-phase periodic solution is parameterized by two pairs of complex conjugated zeroes 
Aj of the polynomial (see [2311201) 



P(A) = J|(A - A,) = A^ - siA^ + s^X^ - S3\ + s^, 
1=1 

Ai = a + i7, X2 = P + iS, A3 = a — 17, A4 = — 16, 
so that the intensity of the wave is given by 

u = |M(x,t)p = z/3 - (z/3 - z/2)sn2(V«^3 - z/i6',m) 



where 



= (7 + 5)2 - 475 sn2( V(a - /3)2 + (7 + 5)2 e, m), 
z/3 - z/2 (Ai - A3)(A2 - A4) 475 



z/3 - (Ai - A4)(A2 - A3) (« - /?)2 + (7 + 5)2^ 

'^i, ^^2, ^3 are the zeroes of the resolvent 7^(z/) of the polynomial -P(A), 



i^i = --(Ai - A2 + A3 - A4)2 = -(a - 



(67) 
(68) 

(69) 
(70) 



z/2 = -^(Ai - A2 - A3 + A4)^ = (7 - 5)^ 
z/3 = -^(Ai + A2 - A3 - A4)' = (7 + Sf, 



(71) 
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and V satisfies the equation 

3 

ve = 2^-n{u), n{u) = Y[{u-u,), z/2<z/<z/3, (72) 

1=1 

9 being the phase of the wave 

4 

e = x-Vt, V = -si = - (73) 

i=l 

In this case there are two fields u and u* , and the corresponding perturbations are equal 

to 

R = -eu, R* = -eu\ (74) 
where £ is a small constant parameter. Easy calculation gives 

dA „ dA OR OA d^R OA ^ , ,2 , 
TT^ + 7)—^ + ^ — ITT + TTT^ = = 2^^^' 

OU OUx ox OUxx OX'^ OU* 

so that we have to average the expression 2eh'g. Now, g{x, t) = X — /i(x, t), where /i satisfies 
the equation 

fie = (75) 
and can be expressed as a function of u (see [TBI 1^ ISO]). 



s g + iV-7^(z/) 
M^) = 4 , ^2<^<^3, ^76) 



Thus, we obtain 



{2evg) = — / v{\i - ^)dx = -(^z/ Ai-- + ^ ^ 



L Jo ' ' 2v ) ^^n{u) 

2e p (Ai - s/4)z/ + g/2 



(77) 



L a/(z^3 - Z/)(z^ - V2){l^ - 1^1) 



where 



^ _ 1 /• d/i _ r« dz/ _ 2K(m) _ 2K(m) 

~ ^ T^TO ~ 7^W^ ~ ~ V(Ai-A4)(A3^' 



2 



The integral in (j77j) reduces to standard elliptic integrals, 

{2sug) = ^ jf {[(2A, - s/2)u^ + q] K(m) + (2Ai - s/2){u^ - z/i)E(m)} . (79) 

-^V ^^3 ~ '^l 

The averaged values (1/(7) and {B/g) are calculated as in the previous cases of the KdV 
equation. As a result, we arrive at the perturbed Whitham equations in the form 

dXi dXi , , 

-^ + f*^ = -^P^, ^ = 1,2,3,4, (80) 
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where Whitham velocities are given by the expression 



-s. + -^=(l-^a.]V, ^ = 1.2.3.4. 



(l/g) dL/dX, V d,L 

which reduces after substitution of (fTSj) to well-known formulae fSE], and functions pi in 
the right hand side of Eq. (jHUj) are equal to 



Pi = 

P2 
P3 

Pa = 



[(2Ai - s/2)ui + q] K(m) + (2Ai - s/2){i^-i - i^i)E{m) 

(Ai - A4)[(Ai - A2)K(m) + (A2 - A3)E(m)] 
[(2A2 - s/2)u^ + q] K(m) + (2A2 - s/2){u3 - u^)E{m) 

(A2 - A3)[(Ai - A2)K(m) - (Ai - A4)E(m)] 
[(2A3 - s/2)u, + q] K(m) + (2A3 - s/2){us - Ui)E{m) 

(A2 - A3) [(A3 - A4)K(m) - (Ai - A4)E(m)] 
[(2A4 - s/2)ui + q] K(m) + (2A4 - s/2)(z/3 - i^^)E{m) 

(Ai - A4)[(A3 - A4)K(m) + (A2 - A3)E(m)] 



^2) 



Equations (jSTj) determine evolution of the Riemann invariants Aj, i = 1,2,3,4, due to com- 
bined action of modulation and perturbation effects. 
For perturbed defocusing NLS equation 



iut + Uxx + 2|'u| u = —ieu (83) 

calculations are actually the same, so we present here only the final results. The intensity 
of periodic wave is given by the formula 

u = |m(x, t)p = z/i + (z/2 - z>i)sn^(A/z/2 - ui 6,m) 

= \iXi -X2-X3 + A4)' + (A2 - Ai)(A4 - A3) sn^ (v/(A3-Ai)(A4-A2)^,m) ^^^^ 

where the zeroes Aj, i = 1,2,3,4, of the fourth degree polynomial -P(A) are real Riemann 
invariants, Ai < A2 < A3 < A4, connected with i/j, z = 1, 2, 3, by the formulae 



z/i = |(Ai - A2 - A3 + A4)^ z/2 = ^(Ai - A2 + A3 - A4)^ 

^^3 = ^(Ai + A2 - A3 - A4)^ 



^5) 



and 



^ ^ i?2-iyi ^ (A2- Ai)(A4-A3) ^ .gg. 
J^s-J^i (A3 - Ai)(A4 - A2)' 



e = x-Vt, V = ~si = -^Xi. (87) 

The wavelength is equal to 



i=l 



^ ^ 2K(m) ^ 2K(m) 

V^^^i V(A4-A2)(A3-Ai)' 
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The perturbed Whitham equations are given by 



where, as usual, 



and 



^ + t'*^ = -^P., ^ = 1,2,3,4, (89) 



(90) 



Vi = 


(1/q) [ V 


^ = 1,2,3,4, 


[(2Ai 


- s/2)u-i + q] K(m) + (2Ai 


-s/2)(z/3-z?i)E(m) 


[(2A2 


(A3-Ai)[(A4-Ai)K(m)- 
- s/2)z73 + q] K(m) + (2A2 


(A4 - A2)E(m)] 
-s/2)(z73-z>i)E(m) 


[(2A3 


(A4-A2)[(A3-A2)K(m)- 
- s/2)z73 + q] K(m) + (2A3 


(A3 - Ai)E(m)] 
-s/2)(z73-^>i)E(m) 


[(2A4 


(A3-Ai)[(A3-A2)K(m)- 
- s/2)z73 + q] K(m) + (2A4 


(A4 - A2)E(m)] 
-s/2)(z73-z>i)E(m) 




(A4-A2)[(A4-Ai)K(m)- 


(A3 - Ai)E(m)] 



(91) 



Pi 

P2 
P3 
P4 

Thus, we see that hnear damping in the NLS equation leads to perturbed Whitham equations 
in simple enough form suitable for applications. 

It is easy to check that more general form of nonlinear damping, 

iWi + Ua^x ± 2|m|\ = -i/(|M|^)M, 

(with a polynomial function fi^u)) as well as action of external potential V{x) slowly de- 
pending on X, 

iut + Uxx ± 2|mP'u = V{x)u, 

can be treated in the same way and arising in this cases integrals can be reduced to standard 
elliptic ones, though the formulae become more cumbersome and they will be presented 
elsewhere together with appropriate application of the Whitham equations. 



4 Conclusion 

The examples considered in the preceding section show that the general formulae ()41|) . p^ 
provide an effective method of derivation of perturbed Whitham equations. They can be 
applied to a number of physically important equations integrable by the 1ST method with 
2x2 matrix or second order scalar spectral problem. For one-phase periodic solutions, 
the perturbation terms in Whitham equations are expressed in the form of integrals which 
can be evaluated with the use of known formulae obtained in derivation of the periodic 
solution under consideration. Thus, the developed here theory permits one to treat effects 
of perturbations in various physical problems related with modulated nonlinear periodic 
waves. One such application has been done in the recent paper |3T] to evolution of modulated 
periodic waves in the Bose-Einstein condensate under slow variation of the nonlinear coupling 
constant. 
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Appendix. Derivation of the identity (1341) 

For simplicity we consider a strictly periodic case when Um{x,t) are periodic functions of x 
with period L. We start from the statement that Eq. is satisfied by the function 

^± = v^exp (^±iV^y/' 1^ (A.l) 

provided g satisfies Eq. (j26j) (see, e.g. "32]). These are Baker- Akhiezer functions which can 
be considered as generalizations of Bloch function solutions of the Schrodinger-like equation 
((7|) with complex "potential" which is a periodic function of x with the period L. Then we 
can introduce formally a "quasi-momentum" p as follows: 

,± ^ v^y^^ ["-dx ^, 

ip = exp(±ipa;), p = / — . (A.2) 

L Jo g 

The "quasi-momentum" p is a functional depending on Um, m = 1, . . . ,N, and their space 
derivatives through the coefficients of the expansion (jHT|) . If we change Um a little, then p 
will change a little too, what defines the functional derivative 



6p ^/^yi'^ ^ 



6Ura L SUm \g 



(A.3) 



where S/6um is defined by (jHSj) . 

On the other hand, this functional derivative can be calculated directly in the follow- 
ing way. We take two functions u'^, u'^ with the same period L and write equations for 
corresponding Bloch functions, 

€. = Au'J^', €. = Aul)r, (A.4) 

which give at once 



Now, from Eq. ()A.1|) we find 

^± = l(^.±2iv^A^^/2)^±, (A.6) 
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so that if V'^ has the property 



then has the same property, 



due to periodicity of g. Hence, for ip' = ip' and ip" = we find from (jA.5jl that 

(e'(^"-^')^ - 1) {€-r^ - ^'~v^r).=o = £ i^i^'J - Au'D) i^'-r^dx. 

Taking the difference — u'^ = 6um small, we obtain 

and for the right hand side of Eq. ()A.9|) we get with the same accuracy 



dA ^ d f dA ^ 



-9 



dUm. dx KdUrr,.,; 



-9 



6umdx, 



(A.7) 



(A.8) 



(A.9) 



so that Eq. ()A.9;i yields 



Sp 

5Um. 



1J-0 



& f dA ' 
9 I , 



du 



(0- 



(A.IO) 



where Am is the order of the highest derivative Um""^ in A. Then comparison of ()A.2|) and 
()A.10|) yields the desired identity 
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